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The ALPHA collaboration has determined the 0(a) im- 
proved Wilson quark action for lattice spacings a < 0.1 fm, in 
the quenched approximation. We extend this result to coarser 
lattices, a < 0.17 fm, and calculate the hadron spectrum on 
them. The large range of lattice spacings obtained by com- 
bining our results with earlier ones on finer lattices, allow 
us to present a convincing demonstration of the efficiency of 
non-perturbative 0(a) improvement. We find that scaling vi- 
olations of the hadron masses studied drop from 30 — 40% for 
the unimproved Wilson action on the coarsest lattice to only 
2 - 3%. 

PACS: 12.38.Gc, 11.15.Ha 



Introduction. To measure standard model parame- 
ters, like CKM matrix elements and quark masses, and 
to find signatures of new physics, accurate knowledge 
of weak matrix elements between hadronic states is re- 
quired. Lattice QCD is the only systematically improv- 
able method of obtaining this information. The high cost 
of lattice QCD simulations has lead to a renewed appre- 
ciation of the fact that progress in this field depends to a 
large extent on the successful use of "improvement" ideas 
(see the proceedings of the last few Lattice Field The- 
ory conferences, e.g. for the last one). The reason is 
the following. To avoid doublers, the Wilson-type quark 
actions most commonly used in simulations must break 
chiral symmetry at some level. On the quantum level at 
least, this violation will generically occur at leading or- 
der in the lattice spacing, 0(a). These errors therefore 
decrease only slowly with the lattice spacing and their ab- 
solute value is large, as experience has shown. To perform 
accurate and reliable continuum extrapolations would re- 
quire the use of very fine lattices, for which simulations 
are very expensive. 

A much better approach Q is to correct the dis- 
cretization errors of a lattice action by adding higher- 
dimensional (irrelevant) operators to the action which 
reproduce the effects of the UV modes omitted on the 
lattice. Trying to do so perturbatively did initially not 
appear to be a significant improvement. It was then re- 
alized that large perturbative corrections arise due to 
lattice-specific "tadpole" graphs, and can be corrected 
by a mean-field type method. Nevertheless, as a resum- 
mation of certain graphs in perturbation theory, this ap- 
proach can basically only reduce quark errors from 0(a) 
to order g^a or g'^a. This is only a logarithmic suppres- 
sion compared to 0(a) and would still require the inclu- 
sion of at least g^a (say) and terms in an honest con- 
tinuum extrapolation of the discretization errors. This 
leads to large errors and potentially unstable fits. 



To eliminate the 0(a) errors of spectral quantities 
there is only one term that has to be added to the Wilson 
QCD action [Q. The gauge action retains the standard 
plaquette form, and the quark action (density) becomes 
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Here and are the standard covariant first, respec- 
tively, second order lattice derivatives. The new a ■ F 
term involves the cr- matrices a^v = — f Iv] and a dis- 
cretization of the field strength F^^ . Inspired by the form 
of its most popular discretization, this term is also known 
as the "clover" term, and the coefficient w as the clover 
coefficient. To eliminate 0(a) errors, uj has to be deter- 
mined as a function of the gauge coupling g. 

A great step forward was recently taken by the AL- 
PHA collaboration , which used the chiral Ward iden- 
tity as an improvement condition to determine the non- 
perturbative value of Lo. This was accornplished in the 
context of the Schrodingcr functional |6|, where one 
imposes fixed boundary conditions on the gauge and 
fermion fields in the time direction, and can then work 
at zero, or at least small, quark masses. The AL- 
PHA collaboration determined improvement coefficients 
for lattice spacings of about a < 0.1 fm (more precisely, 
f3 = &/g'^ > 6.0 in standard notation). 

Since one needs a minimum of three or four reasonably 
separated lattice spacings to perform accurate and reli- 
able continuum extrapolations, this goal will not easily 
be accomplished, even in the "quenched" approximation 
(where quark loops are ignored, and to which the above 
results refer), if only lattices of spacing 0.1 fm and less are 
considered. We will explicitly see this below. We have 
therefore attempted to extend the results of the ALPHA 
collaboration to coarser lattices. 

Chiral Symmetry Restoration at 0(a). Con- 
sider QCD with (at least) two flavors of mass-degenerate 
quarks. The idea [|| for determining the clover coefficient 
is that chiral symmetry will hold only if its Ward iden- 
tity is satisfied as a local operator equation. In Euclidean 
space this means that the PCAC relation between the 
iso- vector axial current and the pseudo-scalar density, 
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should hold for all operators O, boundary conditions, x 
(as long as x is not in the support of O), and also for 
volumes that are not necessarily large in physical units. 
More precisely, it should hold with the same mass m up 
to errors. This will only be the case for the correct 
value of the clover coefficient. 
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Several issues have to be addressed before this idea 
can be implemented in practice. First of all, even though 
here we can ignore the multiplicative renormalization of 
A'' and P^, there is an additive correction to A'' at 0(a), 
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The determination of uj is therefore tied in with that of 
the axial current improvement coefficient ca- Since in 
principle (^) provides infinitely many conditions, this is 
not a fundamental difhculty. How to solve it in practice 
is discussed in 

Note that lo and ca have an 0(a) ambiguity; different 
improvement conditions will give somewhat different val- 
ues for w and ca- Instead of assigning a systematic error 
to Lu and ca one should choose a specific, "reasonable" 
improvement condition — the difference in observables 
from this versus some other choice is guaranteed to ex- 
trapolate away like O(a^) in the continuum limit. 

For various reasons it is preferrable to impose the 
PCAC relation at zero quark mass. Due to zero modes 
this is not possible with periodic boundary conditions; 
the quark propagator would diverge. Another reason to 
abandon periodic boundary conditions is that to be sen- 
sitive to the value of lu it would be highly advantageous 
to have a background field present; it couples directly to 
the clover term. 

The Schrodinger functional provides a natural setting 
to implement these goals. By choosing suitable boundary 
conditions at the "top" {xq = T) and "bottom" (xo — 0) 
of the lattice world, one induces a chromo-electric classi- 
cal background field, and, at least at weak coupling, the 
quark operator has no zero modes at vanishing quark 
mass (the lowest eigenvalue being of order 1 /T) . 

We must now choose a specific improvement condi- 
tion for uj. The idea is that by averaging Eq. (^ over 
spatial volume, each choice of O defines an estimate 
'nio{xo) of the current quark mass. Requiring the dif- 
ference ATO(a;o) = mo^{xQ) — mo2{xo) for two specific 
Oi and O2 to vanish for suitable xq, provides a non- 
perturbative condition to fix u. In practice, one calcu- 
lates all required correlation functions in a Monte Carlo 
simulation for several trial values of lu and finds the zero 
crossing of Am{xo) (more precisely, one should equate it 
to its small, order tree- level value). This determines 
the non-perturbative lu, with some statistical error, for 
the chosen value of the gauge coupling. 

A natural choice of Oi and O2 is provided by boundary 
fields § associated to the lower and upper boundaries 
of the lattice. We will not elaborate on these and other 
choices one makes in the calculation of lu; the details have 
been discussed in the literature [||j7j and the specifics of 
the simulations described here can be found in j|] . 

We have to mention, however, one important point. 
The above simulations at different trial values of lu should 
be performed at a fixed value of the quark mass (de- 
fined by, say, m = ma-^{zo) for suitable zq), preferrably 



zero. It turns out that in the quenched approximation 
this is not possible on coarse lattices: Despite the non- 
periodic boundary conditions one finds in practice that 
for roughly /3 < 6.0 one occasionally hits configurations, 
known as "exceptional configurations" , with an acciden- 
tal (near-)zero mode, leading to a (near-)divergence of 
the quark propagator. (With periodic boundary condi- 
tions configurations with near-zero modes at small quark 
mass exist for any finite /3 in the quenched approxima- 
tion; however, their frequency rapidly decreases at weak 
coupling.) They can be avoided by using a larger quark 
mass, but the question is to what extent this affects the 
value of LU. Fortunately, it turns out that the mass de- 
pendence of LU is extremely weak, so that one can reliably 
determine lu at larger masses. This is illustrated in fig- 
ures |l| and H for coarse lattices (cf. also |^). 

For use of the non-perturbatively improved action in 
later simulations it is advisable to present the results for 
LU in terms of a smooth function of the gauge coupling. 
Combining the results of the ALPHA collaboration |^ 
with our measurements for /9 = 5.7, 5.85, 6.0 and 6.2, we 
find that they can be represented by 
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for P = 6/g^>5.7. This curve incorporates the one-loop 
perturbative result |T^]. It never deviates by more than 
1.0% from the curve presented in Q for /? > 6.0. This is 
illustrated in figure |^, where we used the parameteriza- 
tion of the string tension from to present the clover 
coefficient as a function of lattice spacing. 
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FIG. 1. The non-perturbative clover coefficient as function 
of quark mass and volume for (3 = 5.7. We also show our 
choice of the m = value. 
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FIG. 2. As in figure |l] for = 5.85. 
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3. The measured non-perturbative clover coefficient 



and its parameterization for /3 > 5.7 (solid line). The dashed 
line denotes the curve from The tree-level tadpole es- 
timate from the plaquette is also shown (□). (Using the 
mean-link in Landau gauge gives an estimate closer to the 
non-perturbative determination, cf. 

Hadron Spectrum. To check how small scaling vi- 
olations of spectral quantities are after non-perturbative 
improvement of the action, we have calculated the hadron 
spectrum using Eq. (||) for /3 = 5.7 and 5.85. For a scal- 
ing check it is not necessary to consider light hadrons. 
To avoid the uncertainties of the chiral extrapolation we 
will instead consider hadrons at a pseudo-scalar to vec- 
tor meson mass ratio of mp/mv — 0.7, corresponding 



roughly to the strange quark. This also avoids problems 
with exceptional configurations, which afflict simulations 
at smaller masses on our coarsest lattice. We regard them 
as an essentially technical problem of Wilson-type quarks 
in the quenched approximation (it does not occur for full 
QCD or staggered fermions), orthogonal to the issue of 
improvement. 

Masses were obtained through two-exponential fits 
of correlators from one under- and one over-smeared 
source. We used 400 configurations, statistically en- 
hanced through the use of sources constructed by super- 
imposing different origins with random ^3 phases [ p^ . 
Our results are given in table ^. We also show data from 
other groups on finer lattices, which we interpolated to 
mp/mv = 0.7. Since we can not do correlated fits of 
their data, we multiplied the naive error from interpo- 
lating fits with a factor of 1.5. This gives values close to 
the actually measured errors for neighboring mass values. 
We hope that in the future it will become customary to 
quote hadron masses interpolated to mp/my — 0.7 and 
perhaps a few other benchmark values (like 0.6 and 0.5). 
The results in (l|Jll I were obtained using the parame- 
terization of ui fromg, instead of Eq. (^. We estimate 
that this changes the masses by less than 0.4%, which 
is negligible compared to the current statistical errors. 
The string tensions were taken from our interpolation for- 
mula JlH , which is based on recent precise measurements 
by us and others. We assign these string tensions a 
1% (or smaller) error, that can be added at the end. We 
find that excellent fits to a const -I- ansatz are possible, 



yielding my j y/a = 2.351(20) and m^^ j ^/^ = 3.466(36). 

We have also considered joint fits with data for the 
standard quenched Wilson QCD action {uj = 0). Results 
from different groups for seven couplings in the range 
f3 — 5.7 — 6.3 have been conveniently collected in |l^] 
(errors are treated similarly as above). In a joint fit we 
demand that the ansatze for the improved and standard 
Wilson data intercept at the same point in the continuum 
limit. The results are shown in table || and figure ^. The 
joint fits agree perfectly with fits using only the improved 
action data. For the Wilson data it is necessary to have 
0(a) and O(a^) terms in the ansatz to get a reasonable Q 
in fits where /3 = 5.7 is included. Fitting the Wilson data 
alone yields fits that have either bad Q's or large errors; 
they are also not very stable under leaving out small (or 
large) (3 points. This illustrates how difficult it is to per- 
form reliable continuum extrapolations with the Wilson 
action. Figure ^ also demonstrates that continuum ex- 
trapolations using only lattices with /3 > 6.0 (a^cr < 0.05 
or about a < 0.1 fm) would be quite expensive. 

Conclusions. Figure ^ is impressive proof for the ef- 
fectiveness of non-perturbative 0(a) improvement: The 
scaling violations at /? = 5.7 are reduced from 41% to 
3% for the vector meson mass, and from 33% to 2% for 
the "nucleon" mass. Even more important, the scaling in 
figure H indicates that 0(a) errors really have been elim- 
inated from the improved action to high precision. We 
should remark that without the accurate string tension 
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measurements from it would have been impossible to 
reach this conclusion. 

An analysis of the above data and some toy examples 
shows that it is a factor of 100 or so cheaper to achieve a 
1% (say) error in the hadron masses using the improved 
instead of the standard Wilson action. Since there is 
no fundamental difference in the improvement program 
between quenched and full QCD, we expect very large 
improvements also in more realistic situations like full 
QCD with lighter quark masses. 
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TABLE I. Simulation parameters, string tensions , and 
results for the vector meson and "nucleon" (octet) masses at 
mp/mv ~ 0.7 for the non-perturbatively improved action. 
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FIG. 4. The hadron spectrum from Wilson and improved 
actions at mp/mv ~ 0.7. Also shown are joint fits of both 
data sets (the first vector meson, respectively, nucleon fit from 
table 0). 
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TABLE II. Fit parameters and confidence level Q for 
joint and separate fits of the improved and Wilson hadron 
mass data (at mp/mv = 0.7) to ansatze of the form 
mv / \fa = Vb + Vi a^fa -\- V2 a^a (for the vector meson; simi- 
larly for the nucleon). 
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